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We analyze simplicial quantum gravity in four dimensions using the Regge 
approach. The existence of an entropy dominated phase with small negative 
curvature is investigated in detail. It turns out that observables of the sys- 
tem possess finite expectation values although the Einstein-Hilbert action is 
unbounded. This well-defined phase is found to be stable for a one-parameter 
family of measures. A preliminary study indicates that the influence of the lat- 
tice size on the average curvature is small. We compare our results with those 
obtained by dynamical triangulation and find qualitative correspondence. 
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I. INTRODUCTION 



General relativity relates classical gravitation to the curvature of space-time and therefore 
quantum gravity is essentially a quantum theory of space-time geometry In the search 
for such a theory the sum-over-histories approach has proved to be an elegant and powerful 
tool 0,0,0- It leads directly from the classical action to the quantization of space-time 
using the Feynman path integral. The problems of this approach are well known. Its 
perturbation theory is not renormalizable and a unique prescription for the summation over 
all 4-geo metrics does not exist. 

Moreover, in the case of pure gravity the action is unbounded from below due to rapid 
conformal fluctuations. At first sight this means that the corresponding Euclidean path- 
integral is ill-defined and the quantum theory has no ground state. However, there are 
exceptions such as the Hydrogen atom whose quantum mechanics is well-defined although 
the Euclidean action is unbounded from below. In the case of quantum gravity numerical 
simulations within the Regge calculus indicate that the entropy of the system can suppress 
geometries with large curvature leading to a phase with finite expectation values It 
is the aim of this article to further investigate the conditions for the existence of the well- 
defined phase. We do not give a complete picture of the phase structure and the details of 
the phase transition but refer the reader to appropriate literature 

We use the path integral 

Z = j Dge^^^^e) (1) 

as starting point for the non-perturbative investigation of quantum gravity. The functional 
integral extends over a class of closed 4-geometries g in the Euclidean sector with pure 
gravitational action 

- /^(g) = j d'x^iR - 2A), (2) 

where Lp is the Planck length, R the curvature scalar, and g the determinant of the metric. 
In general the action contains the cosmological constant A. 
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Additional prescriptions are necessary to define the path integral as a limit of systematic 
approximations. A direct route to such approximations is the Regge calculus that uses 
simplicial lattices to discretize general relativity in a coordinate independent and thus very 
elegant way P,p!0[|. One considers a simplicial net being the triangulation of a given 4- 
topology and introduces a metric by assigning a length to each link in the net. Increasing 
systematically the number of simplices and decreasing the link lengths one reaches the 
classical continuum limit if certain conditions are fulfilled |Tl|,|T2|JT^ . 

The Regge approach is a valuable tool to approximate the path integral on a simplicial net 
with fixed incidence matrix and varying link lengths [|6|,^ p!3| , p!^ . A complementary method 
using fixed link lengths and varying incidence matrices is known as dynamical triangulation 
||r5| , p^JT7| . However, the severe difficulties of the functional integral (|l|) are also present on 
the simplicial lattice. 

(i) In general the integration should include a summation over all possible 4-topologies, 
but topologies are not classifiable in 4 dimensions |T^JT9|j20| . 

(ii) In the case of pure gravity the action is unbounded due to rapid conformal fluctuations 
P PT|J2^ . This divergence is also present in the Regge action if one does not restrict or even 
fix all link lengths. The introduction of a cosmological constant term guarantees a finite 
volume of the simplicial lattice, but this is only a necessary and not a sufficient condition 
for finite link lengths and a bounded action. In Sect. 2 we give an explicit example of a 
configuration with finite volume but infinite link lengths and action. One has to explain why 
such configurations do not contribute to the path integral in the well-defined phase and this 
is the aim of the entropy argument of Berg: The integral Z = / dln{l)e~^ is well defined 
if the density of states n{I) vanishes rapidly enough for / —>■ — oo. Previous computations 
indicate that the entropy of the system can indeed compensate the unbounded action leading 
to a phase with small and negative average curvature |^,|^ and one has assumed that a finite 
volume will always lead to finite link lengths [^]. 

To give, at least numerically, further evidence for this assumption we investigate in Sect. 
3 the well-defined phase of the Regge approach with three different methods. Furthermore, 
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we compare our results with those obtained by dynamical triangulation and find qualitative 
agreement. 

(iii) A unique definition of the measure does not exist, since different physical ar- 
guments prefer different classes of measures to be outhned below [p|, p3| , ^p5| , ^ . It is an 
open question in what manner various choices of the measure in the path integral affect 
the well-defined phase although one expects universality for a certain class 0,1^]. Hence, 
we use a one-parameter family of measures and investigate in Sect. 4 the stability of the 
entropy dominated region against variations of this parameter. Again, the comparison with 
dynamical triangulation shows a remarkable coincidence of both approaches. 

(iv) The discretization of general relativity destroys the diffeomorphism group and it is 
not known if it can be recovered in the continuum limit |lT0| , p7| . Until now the existence 
of a reasonable continuum limit has not been demonstrated for simplicial quantum gravity. 
Therefore, the meaning of our numerical results for the continuum remains unclear. However, 
it is possible to study the influence of flnite-volume effects on the computed expectation 
values and it seems that the entropy dominated phase survives increasing lattice volumes 
as demonstrated in Sect. 4. This might be an important result for computations within the 
'fundamental lenght scale' scenario proposed recently by Berg et al. |p8|| . 

The article is organized in the following way: The concept of the Regge calculus and 
the simplicial path integral are discussed in Sect. 2. In Sect. 3 we investigate the entropy 
dominated phase by checking various observables for well-deflned expectation values. The 
stability of this phase against a variation of the measure and flnite-volume effects are studied 
in Sect. 4. The conclusion follows in Sect. 5. 



II. SIMPLICIAL QUANTUM GRAVITY 

The triangulation of a manifold Ai is determined if a simplicial manifold /C homeomorphic 
to Ai is given |]TD|J2^. A regular triangulation of the 4-torus used in the following is 
provided by an algorithm that divides into hypercubes and then each hypercube into 
24 4-simplices For the total numbers of d-simplices one has A^o = n'^, n > 3, and 
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Ni = 15No, N2 = 50No, N3 = 6ON0 and A^4 = 24A^o- A triangulation of with a smaller 
number of simplices is described in [^]. In case of the hypercubic triangulation of T*^ the 



coordination numbers Nn/Nm are close to those of a random triangulation as proposed by 
Christ et al. 

A simplicial net obtains a metric structure if one assigns a length to each link of the 
lattice. In the following we denote the squared length of link I hy qi E R}. Since we consider 
the Euclidean sector, we require > and have to check the Euclidean triangle inequalities 
for all simplices P,|9|,p!0[] . In other words, for every 4-simplex s in the lattice we demand 
that it can be constructed in W^. After fulfilling the Euclidean triangle inequalities for every 
4-simplex one has a consistent configuration {qi} that determines the simplicial lattice. 

One can calculate several quantities of a simplicial lattice such as the area At of the 
triangle t = 1, A^2 and the 4- volume Vs of the 4-simplex s = 1, A^4. Another important 
quantity is the deficit angle 6t associated with each triangle t. It is defined as 

6t = 2n-Y,e,,t, (3) 

where the sum runs over all 4-simplices s sharing the triangle t. The interior angle Og^t is 
the dihedral angle of the 4-simplex s with the basis triangle t . 



From the viewpoint of differential geometry a simplicial lattice is a singular 4-geometry 
with the curvature concentrated on the triangles and we consider the simplicial lattice as an 
approximation to a smooth Riemannian manifold |]5|,|6|j9| JlO|JTT|JT^Jl3|l . The connection of the 



deficit angle with the curvature tensor is exhibited by parallel transporting a vector along a 
closed loop around a triangle t P,P,|5^ . 

In consideration of the proportionality between the curvature tensor and the deficit angle 
and guided by the dimension of the action Regge proposed [^] 

/ d'^x^R ^ 2j2At6t. (4) 



An alternative derivation and an action with the same continuum limit can be found in 32 



One can show that the Einstein-Hilbert action of a smooth 4-geometry converging to- 



wards a simplicial lattice approaches the Regge action [0, but for our purposes it is more 



important to study the opposite limit, approximating a given smooth 4-geometry by increas- 
ing the number of simphces of the lattice and decreasing the link lengths systematically. A 
necessary prerequisite in this discussion is the definition of the fatness (ps of a 4-simplex s. 
It is a relation between the maximal link length and the 4-volume given by 

0s = ^V^, (5) 

max;e^(g,4) 

where the constant is set to C = 24 in the following. (Cheeger et al. use a slightly different 
but equivalent definition |[13|.) If the lattice skeleton approaches a smooth geometry such 



that 0s > / = const > one recovers general relativity in the continuum limit []T2|,|TB|. 



Notice that the Regge calculus is the only known approximation scheme of general relativ- 
ity incorporating the Bianchi identities • Since these relations connect the dynamical 
and constraint equations of general relativity they are essential for any attempt to quantize 
gravitation. In the following we try to assign a concrete meaning to the path integral (|l|) in 
the language of the Regge calculus. 

(i) The integration symbol / stands for a summation over all 4-geometries {Ai, g) without 
boundary. This can be replaced by a sum over distinct simplicial lattices {IC,{qi}). One 
usually restricts the class of incidence matrices to include only simplicial manifolds /C with 
a fixed topology. Further simplifications lead to two somehow complementary methods. 

Dynamical triangulation: One assigns the same length to each link, qi = IV/, and 
performs a summation over all possible triangulations /C of a given topology 

Z = J2e-'^^^^. (6) 

K. 

The main motivation for this approach is its success in two dimensions . In four 



dimensions the action simplifies for equal link lengths to —Ie = C2N2 — C4N4 with coupling 
parameters C2, C4. It is bounded from below (and above) for a finite number of simplices and, 
moreover, since the inequality N2 < ION4 holds for every simplicial lattice one has Ie > 
for C4 > 10c2. This can be fulfilled by making the cosmological constant large enough, in 
contrary to the continuum theory and the Regge approach. Since we try to investigate the 
possible existence of an entropy dominated phase for unbounded gravitational action, we 
use the 
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Regge approach: The triangulation is kept fixed and the path integral Z reduces to a 
summation over different configurations {qi}- Only the link lengths are allowed to vary but 
not the incidence matrices [^|8|, p!0| , p!^ , 



Z = f Dge-^^(^^'». (7) 



Important motivations for this approach are studies of the weak field limit which show the 
correspondence with linearized quantum gravity Comparisons of Regge calculus and 
mini-superspace approximations have also been performed with rather promising results 

MM- 

(ii) The gravitational action /^(g) is unbounded in the continuum due to conformal 
fluctuations and this unpleasant feature is also present in the Regge approach yet for a 
finite number of simplices • If some of the 4-simplices collapse the corresponding areas 
can grow, — oo, although their 4- volumes stay finite. To give an example consider a 
configuration {qi} where first all link lengths are equal to a, i.e. qi = VZ. Then choose 
one vertex vq of the lattice arbitrarily and assign to each link of this vertex a length L > a. 
It is easy to check that the Euclidean triangle inequalities are all satisfied in this case. Now 
consider a ^ and L — > oo so that a^L = const. In this case the 4-volume of the 4-simplices 
with vq as a vertex vanishes like a^L and since the 4-volume of the other simplices disappears 
as the total 4-volume of the lattice approaches zero, J2s K ~^ 0- contrast, the area of 
triangles with vq as a vertex diverges as aL and therefore the action grows without limit, 
J2tASt oo. 

As mentioned in the introduction, it seems that such configurations do not contribute 
to the path integral in the well-defined phase and we investigate this possibility in 

the following section. Of course, the described divergence cannot be present for dynamical 
triangulation with all link lengths equal. 

(iii) A unique definition of the gravitational measure Dg does not exist since it is not 
clear which quantities have to be identified with the 'true' physical degrees of freedom 
§JT|j2|j2|j2D. Within the Regge approach it seems natural to translate the gravitational 



measure into a product over the links of the simplicial lattice 



Dq = '[[Kl)dqiJ^{qi,..,qNi), 



(8) 



where the function ensures that the computations are performed in the Euchdean sector 
being equal 1 if the generahzed triangle inequalities are fulfilled and otherwise. The weight 
function is a (perhaps complicated and non-local) function of {qi} associated with each 
link / [pSj^ . The simplest choice is the uniform measure = 1 VZ and one has argued 
that it corresponds to the DeWitt measure |0,[lO|. However, considering the results of Misner 
and Faddeev-Popov it seems natural to choose so that a scale-invariant measure results 



Relying on this background the presented numerical studies have been performed with 



where the parameter cr > determines the behavior of the measure under rescaling. The 
case (7 = 1 corresponds to the uniform measure while for o" = a scale-invariant measure in 
its simplest form results. 

Metropolis simulations within the Regge approach proceed in the following way [p|,|6i|8|,p!4| . 
Starting with a Euclidean configuration {qi}Q a Markov chain of configurations {qi}T is 
generated to approximate the canonical ensemble. Proposing a new configuration one has to 
check first the Euclidean triangle inequalities. If the new configuration is not in the Euclidean 
sector it is rejected immediately. Otherwise one has to check the Metropolis condition to 
fulfill the requirement of detailed balance. Given a Euclidean configuration one changes 
usually only one link length at once which simplifies the test of the triangle inequalities. 
Scanning one time through the whole lattice is called a sweep. After thermalization of 
the system the averages over the generated configurations are taken to approximate the 
expectation values. 

A priori one expects that simulations of simplicial quantum gravity with unbounded ac- 
tion might not find an equilibrium. Indeed, one has to face this problem for computations on 
non- regular triangulations P^ , |57|] . To investigate systematically if the entropy of the system 
can stabilize the unbounded action we use the following method: Guided by the analysis 



/^(O = qi 



cr-l 



(9) 
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of Cheeger et al. [T^ we set a lower limit for the fatness of each 4-simplex, 0s > / > 0, 
restricting the configuration space. In this way we introduce a scale-invariant cutoff allow- 
ing the system to reach an equilibrium after a finite number of sweeps. As long as / > 
the gravitational action is bounded for a finite number of 4-simplices. Decreasing stepwise 
the value of / we investigate the limit / ^ and conclude from the convergence of this 
process if a well-defined phase does exist or not. Since decreasing / enlarges the significant 
configuration space one has to increase the number of sweeps to ensure equilibrium. In the 



actual computations the number of sweeps varies from 5k (/ = 10 ^) up to 20k (/ = 10" 



III. ENTROPY VERSUS ACTION 

In the following we employ the action in the form 

-iE = pJ2^tSt-\Y.Vs, (10) 

t s 

with a bare cosmological constant A fixing the scale and we restrict our attention to the 
simplicial measure (^) with the parameter cr. Considering the behavior of the path integral 
Z{P, A) under rescaling qi rqi one obtains 

after setting r = j3/\. Differentiating InZ with respect to {3 and A gives 

-(3{Y, AA) + 2A(^ K) = cjN, (12) 

t s 

in the well-defined phase 0. We therefore set the parameter A equal to a so that {Vs) takes 
the value Ni/2N^ for a > at = 0. For /? > 0, (K) decreases slightly as well as {qi). To 
compensate the changes of the lattice size as a function of (3 and a we use the dimensionless 
quantity 

^^ = f(^.) (13) 

to compare the results. The factor | respects the original form of the Regge action (^) and 
since [3 corresponds to 2L~p^ the dimensionless observable £^ expresses the average squared 



lattice spacing in units of the bare Planck length. Different to conventional lattice field- 
theory the simphcial lattice itself is the quantum object, i.e. the lattice spacing is an 
observable rather than a parameter and different definitions are possible §]. 

Another important observable of simplicial quantum gravity is the scale invariant average 
curvature of a lattice configuration |PJ^ 

The expectation value (R) can be understood as effective cosmological constant measured 
in lattice units. This can be seen from the classical field equations deriving A = and 
obtaining 

A = l ^ff^V^R ^ 1^ (^5) 
4 Jd^x^ 2 ' ^ ^ 

with A the classical cosmological constant, R the curvature scalar and g the determinant of 
the metric tensor. 

In Fig. 1 the behavior of the expectation value (R) as a function of i is depicted. The 
most interesting region (R) < is displayed in more detail in the lower plot. The value of 
the lower limit / for the fatness decreases as follows: / = 10~", n = 3,4,5,6,oo. For every 
/ we increase f3 stepwise leading to increasing i. Since the significant configuration space 
becomes larger with smaller values of / we performed 20k sweeps for each data point for 
/ < 10"^ and 5k sweeps otherwise. Averages are taken over the last 10k and 2.5k sweeps, 
respectively, so that the statistical error is smaller than the symbol size. This computation 
was performed with the uniform measure (cr = 1) on the regular hypercubic triangulation of 
the 4-torus with 4"^ vertices. We will see later that the well-defined phase is stable against 
changes of the parameter a in a certain range. 

For all values of the cutoff / we find a well-defined phase with small and negative curva- 
ture. Increasing i one suddenly enters a region of large positive curvature. For / < 10^^ it 
is difficult to reach a stable equilibrium across the transition point due to a trend towards 
non-reproducible states. This coupling regime is therefore called the ill-defined phase. 
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In every case it seems that {R) < as long as £ < ic, with ic = 0.45(1) for sufficiently 
small /. Furthermore, we observe for / ^ a convergence of the function {R){i) at every 
i < ic'- As long as £ < 0.3 one sees practically no difference in the average curvature (R) for 
all values of the cutoff /. This is a strong indication that configurations with small fatness 
(ps give almost no contribution in this region. Increasing i further one begins to realize 
differences in (R) for larger /. However, for / < 10^^ there are (nearly) no differences to 
see over the whole well-defined phase < i < ic- This convergence for / indicates 
that the entropy of the system is indeed able to compensate the unbounded gravitational 
action. Configurations with large average curvature occur with small probability and give 
(practically) no contribution to the path integral in the well-defined phase. Therefore, it 
might be possible to perform computations even without any restriction (/ = 0) and obtain 
reasonable expectation values. One should be aware that the scale of the curvature radius 
should be larger than the average lattice spacing and smaller than the size of the system, 
{q,)<{q,)\{R)\-'<N',/'{qi) 

To further analyze the probability of 'crumpled' configurations with collapsing simplices 
we display the behavior of the average fatness {(ps)- Fig. 2 shows {(ps) as a function of i for 
decreasing cutoff /. The parameters of the system and the configurations sampled for the 
statistics are the same as before. It turns out that (0s) is much larger than the lower bound 
/ in the well-defined region. The convergence of the function {(j)s){i) for decreasing / — is 
seen clearly in the well-defined phase. With the transition to large values of {R) the average 
fatness decreases suddenly as expected. This indicates the collapse of 4-simplices and the 
transition to a 'crumpled' lattice. 

Returning to the well-defined phase let us investigate the behavior of the expectation 
value (6^). The following argument will supply us with a further indication for the stability 
of the well-defined phase. Imagine that one incorporates an additional term —a J2t 
action. Differentiation of In Z then gives 

dlnZ ^ ^, d\nZ /v^ r2\ 
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d{AA) 



(16) 



8(3 



da 



All computations are performed at a = but we assume that pushing a slightly to a positive 
value would shift {At5t) towards zero since the additional term prefers flat configurations 
and supresses the curved. (Recall the definition of R to see that i? < corresponds to 
^t^t^t < 0). This means that {At6t) will increase for negative values and decrease for 
positive. We thus have -^{At5t) > and therefore -§p{5t) < in the well-defined phase with 
negative curvature. Since £ is proportional to the fact that {At5t) is negative induces 
that (6^) decreases in this phase with £. As seen in Fig. 3 {6f) indeed decreases as expected 
with the minimum just around the transition point. This is an indication that the average 
curvature is really negative in this phase and at least shows the self consistency of the 
numerical results. 

It is practically impossible to determine a systematic error for Metropolis simulations. 
One cannot exclude for example the possibility that the well-defined phase is only a numerical 
artefact due to metastable states preventing the system from reaching configurations with 
large positive curvature. To check this we performed computations with inhomogeneous 
start configurations from the ill-defined phase having large positive average curvature. 

Without a restriction of the fatness (/ = 0) and a non-zero coupling (3 the well-defined 
phase i < ic was reached after a large number of sweeps. The average curvature returned to 
small and negative values of R although the unbounded gravitational action prefers a large 
positive curvature. The history of such a run is depicted in the upper plot of Fig. 4 and 
compared with the history for a start with a homogeneous configuration. The independence 
of (R) from the start configuration is the strongest indication that entropy is indeed able to 
compensate the unbounded action in the well-defined phase. 

Near the transition point very long runs are necessary to reach the equilibrium for very 
inhomogeneous start configurations. Limited resources therefore do not allow to perform the 
above check for Rstart oo and i —>■ ic- Thus, we cannot exclude the possibility that near 
£c metastable states occur if the phase transition is first order as reported in [0. However, 
applying a (small) lower limit for the fatness (/ = 10~^) we found no influence of the start 
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configuration for £ < 0.45 witli ratlier large Rstart- Of course, the curvature of the start 
configuration is hmited after applying a lower limit on the fatness. The lower plot of Fig. 4 
shows such a history reaching the equilibrium after 50k sweeps. 

The two complementary approaches to simplicial quantum gravity are compared in Fig. 
5 by plotting where we depict our results within the Regge approach and those obtained by 
Briigmann within dynamical triangulation [^]. The lower curve gives {R) as a function of 
£ for / = 10~^ and uniform measure taken from Fig. 1. The upper curve shows the results 
of Metropolis simulations performed by Briigmann dynamically triangulating the 4-sphere 



with = 16k IQ. We depict (R) as a function of i and since = 1 V/ we have 

^ = ' ^^^^ 

where a,v,6 denote the area, 4- volume and interior angle associated with equilateral sim- 
phces of link lengths equal to 1 and i = \J\^- 

For a given number of 4-simplices the number of triangles in the lattice is limited. Since 
N2 < 10A^4 there is an upper limit for the average curvature, R < (a/f )(207r — 10^), different 
to the Regge approach. However, it seems that both dynamical triangulation and the Regge 
approach lead to similar results at least qualitatively. (One has to remember that the 
underlying topology is different.) The average curvature begins for small C. with small positive 
values and undergoes a transition to large positive curvatures. The transition point seems to 
be located close to that of the Regge approach, but it is unclear if this is pure coincidence. 



IV. INFLUENCE OF THE MEASURE AND FINITE- VOLUME EFFECTS 

In this section we address first the question of the measure in the well-defined phase of 
simplicial quantum gravity and investigate then the influence of finite-volume effects. 

Fig. 6 shows the expectation value {R) as a function of ^ for measures with different 
values of a increasing stepwise from to 1.5 including the scale- invariant and the uniform 
measure. For the scale-invariant measure (a = A = 0) we use directly the constraint 
X)s Vs = const rescaling the lattice at every sweep as proposed by Berg |^. The computations 
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have been performed on the regular hypercubic triangulation of the 4-torus with 4^ vertices. 
To facihtate numerical simulations a lower bound / = 10~^ has been applied for the fatness 
(ps of each 4-simplex. The results of Sect. 3 suggest that such a restriction of allowed 
configurations practically does not affect the behavior of the system in the well-defined phase. 
Averages are taken over at least 10k sweeps after thermalization so that the statistical error 
is smaller than the symbol size. 

For 0" < 1 and < £ < 0.3 the value of (i?) seems to be practically independent of 
cr. Such a stability against a variation of the measure was reported first by Hamber in the 
context of theory [S. A direct comparison of scale- invariant and uniform measure for 



pure gravity with a cutoff for the link lengths is described in Ref. Near the transition 
point, £ ~ 0.4, the influence of a becomes more pronounced. 

The result for a = 1.5 differs over the whole range of ^ from those obtained for cr < 1. 
At £ = the value of {R) is significantly larger than for a < 1. With increasing ^ we find 
a stable equilibrium even for {R) > 0. This suggests that the considered measures fall into 
two qualitatively different classes: While we see almost no dependence of the expectation 
values in the interval < cr < 1 the influence of the measure becomes visible for cr > 1. The 
situation for £ = is further illustrated in the lower plot of Fig. 6 showing (R) as a function 
of a. Setting cr = 2 we found no equilibrium even at (3 = 0. 

To gain more information about the lattice geometry the behavior of the areas and deficit 
angles is examined separately. Fig. 7 displays the scale-invariant quantity {At) / {qi) as a 
function of £ for the different values of cr. For < i < 0.4 this ratio stays almost constant 
somewhat below the value ^ ~ 0.433 of equilateral triangles. The value for a = 1.5 lies 
significantly below the other data points. Across the transition at ^ ~ 0.45 the ratio {At) / {qi) 
decreases, indicating a distortion of the triangles. 

The expectation value of the average deficit angle {6t) is depicted as a function of i in Fig. 
8. Surprisingly {6t) stays negative for all values of i even after the transition to large positive 
curvature. In the well-defined phase the absolute value of {6t) is rather small compared to vr. 
One can understand this behavior by examining the hypercubic triangulation of the 4-torus 
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that contains two different types of triangles. The number of 4-simphces sharing a triangle of 
the first type is 6 while it is 4 for the second type. As noticed first by Berg the contributions 
of these two types almost cancel each other leaving a small negative average deficit angle 
i- 

The dependence of the simplicial path integral on the measure has been studied also 
within the framework of dynamical triangulation. To investigate the influence of different 
types of the gravitational measure Briigmann considered an additional term nJ2v^^{o{v)) 
in the action where o{v) is the 'order' of the vertex v, i.e. the number of 4-simplices that 
contain this vertex |^^. Although the correspondence with the continuum is not entirely 
clear, it is plausible that this term reproduces a measure of the form fi'"^^ ^i<j ^Qij with 
n = —5 the scale-invariant and n = the uniform measure of DeWitt. Results of these 
computations are depicted in Fig. 9 where again the expectation value {R) is plotted versus 
i for different values of n. Since Briigmann considers negative couplings /3 < we extend i 
to negative values. 

Besides a shift in i for dynamical triangulation, a comparison of Figs. 6 and 9 shows 
a remarkable qualitative coincidence of the results with the choice n = for dynamical 
triangulation corresponding to the case cr = 1 of the Regge approach. Again the influence 
of the measure seems to be small between the scale-invariant and the uniform measure, but 
turns out to be significant for n = 5. 

As for every lattice field theory it is important to study finite- volume effects for simplicial 
quantum gravity . To investigate the stability of the well-defined phase for larger lattices 
we increased the number of vertices from 3^ to 8^. The results are displayed in Fig. 10 and 
show {R) as a function of i computed on the hypercubic triangulation of the 4-torus with 
the uniform measure {a = 1) and a lower limit / = 10~^ for the fatness. For the smaller 
lattices with Nq = 3^, 4^, 6^ we performed at least 10k sweeps after thermalization. Statistics 
is rather poor in the case of Nq = 8^ where the number of sweeps in equilibrium is only 2.5k. 
Therefore, one has to consider the three data points for the 8^-triangulation as preliminary. 

First results indicate that the critical value of the coupling (3 and therefore ic slightly 
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decreases for increasing A^o- But it seems that increasing the number of vertices No has 
only a weak influence on the computed expectation values of one-point functions and the 
well-defined phase survives on larger lattices. 



V. CONCLUSION 

The Regge calculus provides a direct route to systematic approximations of the quantum- 
gravity path-integral. It allows numerical studies of non-perturbative quantum gravity and 
has led to astonishing results. 

The most important finding is the existence of a well-defined phase. We studied this 
phase by applying a lower bound on the fatness of each simplex and decreasing it system- 
atically. The convergence of this process supports the 'entropy hypothesis'. Although the 
gravitational action is unbounded due to conformal fluctuations the entropy of the system 
seems to stabilize the expectation values. In a certain range of the gravitational coupling a 
phase with small negative average curvature occurs within the Regge approach. The pre- 
dicted behavior of the squared deficit angle and the independence of the simulation from 
the start configuration further support the stability of the entropy dominated phase. 

The well-defined phase turns out to be stable against variations of the measure and an 
increase of the lattice size. The influence of the measure has been studied using a one- 
parameter family of simple local functions. However, a wider class of different measures 
should be investigated including also non-local types and non-regular triangulations. The 
dependence of the results on the lattice size has been studied to some extent. The fact 
that the entropy dominated phase survives on larger lattices is encouraging for studies of a 



'fundamental length scale' scenario p8 . 



ACKNOWLEDGMENT 

We are grateful to Bernd Briigmann for transmitting the data of his simulations enabling 
the comparison of Regge approach and dynamical triangulation. This work was supported 



16 



in part by " Fonds zur Forderung der wissenschaftlichen Forschung" under contract P9522- 
PHY. 



[1] Ch. Misner, K. Thorne, and J. Wheeler, Gravitation (Freeman, San Francisco, 1973). 
[2] Ch. Misner, Rev. Mod. Phys. 29, 497 (1957). 

[3] S. Hawking, in General Relativity, An Einstein Centenary Survey, edited by S. Hawking and W. 
Israel (Cambridge University Press, Cambridge, 1979). 

[4] J. Hartle and S. Hawking, Phys. Rev. D28, 2960 (1983). 

[5] H. Hamber and R. Wilhams, Nucl. Phys. B248, 392 (1984); Phys. Lett. 157B, 368 (1985); Nucl. 
Phys. B269, 712 (1986). 

[6] H. Hamber, in Critical Phenomena, Random Systems, Gauge Theories, Proceedings of the Les 
Houches Summer School, Session XLHI, edited by K. Osterwalder and R. Stora, Vol.1 (North 
Holland, Amsterdam, 1986). 

[7] H. Hamber, Phys. Rev. D45, 507 (1992); Nucl. Phys. B400, 347 (1993); Nucl. Phys. B (Proc. 
Suppl.) 30, 751 (1993); H. Hamber and R. Williams, Phys. Rev. D47, 510 (1993). 

[8] B. Berg, Phys. Rev. Lett. 55, 904 (1985); Phys. Lett. B176, 39 (1986). 

[9] T. Regge, Nuovo Cimento 19, 558 (1961). 

[10] J. Hartle, J. Math. Phys. 26, 804 (1985); J. Math. Phys. 27, 287 (1986); J. Math. Phys. 30, 452 
(1989). 

[11] R. Friedberg and T. D. Lee, Nucl. Phys. B242, 145 (1984). 

[12] G. Feinberg, R. Friedberg, T. D. Lee, and H. Ren, Nucl. Phys. B245, 343 (1984). 

[13] J. Cheeger, W. Miiller, and R. Schrader, Comm. Math. Phys. 92, 405 (1984); Lect. Notes in 
Phys. 160 (Springer, Berlin, 1982). 

17 



[14] W. Beirl, E. Gerstenmayer, and H. Markum, Phys. Rev. Lett. 69, 713 (1992); Nucl. Phys. B 
(Proc. Suppl.) 26, 575 (1992); W. Beirl, E. Gerstenmayer, H. Markum, and J. Riedler, Nucl. 
Phys. B (Proc. Suppl.) 30, 764 (1993); W. Beirl, E. Gerstenmayer, and H. Markum, in Quantum 
Physics and the Universe, Proceedings of the International Symposium, edited by M. Namiki, 
Y. Aizawa, K. Maeda, and I. Ohba (Pergamon, Oxford, 1993). 

[15] J. Ambj0rn, B. Durhuus, and J. Prohlich, Nucl. Phys. B257, 433 (1985); V. Kazakov, I. Kostov, 
and A. Migdal, Phys. Lett. 157B, 295 (1985); M. Agishtein and A. Migdal, Int. J. Mod. Phys. 
CI, 165 (1990). 

[16] J. Ambj0rn, Nucl. Phys. B (Proc. Suppl.) 25A, 8 (1992); J. Ambj0rn and S. Varsted, Nucl. 
Phys. B (Proc. Suppl.) 25A, 25 (1992). 

[17] J. Ambj0rn and J. Jurkiewicz, Phys. Lett. B278, 42 (1992); M. Agishtein and A. Migdal, Mod. 
Phys. Lett. 7, 1039 (1992). 

[18] R. Geroch and J. Hartle, Found. Phys. 16, 533 (1985). 

[19] J. Hartle, in Proceedings of the Fourth Marcel Grossmann Meeting on General Relativity, edited 
by R. RufHni (Elsevier, Amsterdam, 1986). 

[20] A. Markov, Proc. Int. Congress Math., 300 (1958). 

[21] P. Mazur and E. Mottola, Nucl. Phys. B341, 187 (1990). 

[22] S. Giddings, Int. J. Mod. Phys. A5, 3811 (1990). 

[23] H. Leutwyler, Phys. Rev. 134, B1155 (1964); E. Pradkin and G. Vilkovisky, Phys. Rev. D8, 4241 
(1973). 

[24] B. DeWitt, Phys. Rev. 160, 1113 (1967). 

[25] L. Faddeev and V. Popov, Usp. Fiz. Nauk. 109, 427 (1974); Sov. Phys. Usp. 16, 777 (1974). 
[26] P. Menotti, Nucl. Phys. B (Proc. Suppl.) 17, 29 (1990). 

[27] M. Rocek and R. Wilhams, Phys. Lett. 104B, 31 (1981); Z. Phys. C21, 371 (1984). 



18 



[28] B. Berg, B. Krishnan, and M. Katoot, Nucl. Phys. B404, 359 (1993). 

[29] C. Nash and S. Sen, Topology and Geometry for Physicists (Academic Press, London, 1983). 

[30] W. Kiihnel and G. Lassmann, Discrete Comput. Geom. 3, 169 (1988). 

[31] N. Christ, R. Priedberg, and T. D. Lee, Nucl. Phys. B202, 89 (1982). 

[32] M. Caselle, A. d'Adda, and L. Magnea, Phys. Lett. B232, 457 (1989). 

[33] J. Louko and P. Tuckey, Class. Quantum Grav. 9, 41 (1992). 

[34] B. Briigmann, Nucl. Phys. B (Proc. Suppl.) 30, 760 (1992). 

[35] M. Bander, Phys. Rev. Lett. 57, 1825 (1986). 

[36] A. Jevicki and M. Ninomiya, Phys. Rev. D33, 1634 (1986). 

[37] W. Beirl, Ph. D. Thesis, Technical University of Vienna, 1992 (unpublished). 

[38] D. Gross, Phys. Lett. 138B, 185 (1984); A. Billoire, D. Gross, and E. Marinari, Phys. Lett. 139B, 
75 (1984). 



19 




-0.2 0.2 0.4 0.6 



0.8 

It 



1.2 1.4 1.6 1. 



{R) 




Fig. 1. Expectation value {R) within the Regge approach as a function oii for decreasing lower 
limit / of the fatness of each 4-simplex. Computations were performed on a regular triangulation of 
the 4-torus with 4^ vertices using the uniform measure ((7 = 1). The upper plot gives an overview 
of all results available whereas the lower plot magnifies the well-defined phase. The existence of a 
well-defined phase and the convergence of {R) for / — >^ are clearly seen. The transition point to 
positive curvature is shifted from ic ~ 0.8 for / = 10~^ down to ^c ~ 0-45 for / < 10~^. 
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Fig. 2. Average fatness {(pg) versus £ for decreasing cutoff / for the fatness of each 4-simplex. 
The value of {<ps) stays much larger than / in the well-defined phase and convergence for / — >^ 
is seen. With the transition to large curvatures {(ps) decreases suddenly indicating the collapse of 
4-simplices into degenerate configurations. Configurations and statistics are the same as in Fig. 
1. 
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Fig. 3. Average squared deficit angle ) as function of I for decreasing cutoff /. A non- 
monotonic behavior of {Sf) is seen for increasing I. The fact that the average squared deficit angle 
decreases with increasing £ is a strong indication for the stability of the well-defined phase. Near 
the transition point the minimum deviation from flat space is reached and a further increase of 
/3 and thus £ leads to the ill-defined region with fast increasing squared deficit angles (see text). 
Configurations and statistics are the same as in Fig. 1. 
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Fig. 4. (In) dependence of the numerical results in the well-defined phase from the start config- 
uration. The upper plot shows the histories of R for a start with homogeneous and inhomogeneous 
configurations, respectively. Both simulations converge to the same equilibrium in the well-defined 
phase for non-zero {3 {I k. 0.29). The computer time is denoted by r with one unit representing 
10 sweeps. Although the unbounded action (/ = 0) favors a large positive curvature the entropy 
moves R to negative values. With a small bound for the fatness (/ = 10"^) independence from the 
start is seen even starting with very inhomogeneous configurations as depicted in the lower plot. 



23 



{R) 



40 
30 
20 
10 

-10 



1 1 1 1 1 
-/\ Dyn- Triangulation 
O Regge Approach ^ 

A 

A 
A 


^^^^ 
ill-defined 


well-defined o ■ O ' '* 
1 1 1 1 1 


1 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



e 



Fig. 5. Expectation value (R) versus i for the two complementary approaches to simpli- 
cial quantum gravity. The upper curve shows results from dynamical triangulation obtained by 



Briigmann approximating the 4-sphere by 16k 4-simplices |34|. The lower curve was obtained 
within the Regge approach on a regular triangulated 4-torus with 4^ vertices. The behavior of 
both functions is rather similar although the underlying topologies differ. A phase of small {R) 
turns into a region with large curvature for increasing i. The black dot indicates a value close to 
the transition point ic for the Regge approach. 
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Fig. 6. Expectation value {R) as a function of £ within the Regge approach for different types of 
the measure parametrized by o" > 0. Computations have been performed on a regular tringulated 
4-torus with 4^ vertices applying a cutoff / = 10~^ for faster convergence. The upper plot shows 
that the behavior of {R) in the region of small i is almost independent of a as long as a < 1. 
This stability ceases for cr > 1 as depicted in some detail in the lower plot where (R) is given for 
increasing a in the case of pure entropy, i = P = 0. 
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Fig. 7. Ratio {At)/{qi) as a function of ^ for different measures parametrized by a > 0. In 
the well-defined phase this scale-invariant quantity stays almost constant and somewhat below the 
maximum value ^ 0.433 for equilateral triangles. For the scale-invariant measure {a = 0) the 
upper curve results and the expectation values decrease slightly with increasing a. Configurations 
and statistics are the same as in Fig. 6. 
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Fig. 8. Expectation value as a function of I for different measure parameters a. It is 

remarkable that ((5t) stays always negative even after the transition to positive curvature. The 

curves lie close together for < cr < 1 while the behavior differs significantly for a = 1.5. The 
configurations are those of Fig. 6. 
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Fig. 9. Investigations of different measures with dynamical triangulation of the 4-sphere per- 



formed by Briigmann [34|. An additional term n^^ln(ot,) in the action mimics different proposed 
measures. Besides a shift in £ the picture has a striking similarity to Fig. 6 if one identifies a = 
with n = —5 (scale- invariant measure) and a = 1 with n = (uniform measure). Notice the small 
influence of n in the range —5 < n < +1 and the exceptional behavior of n = +5. 
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Fig. 10. Influence of the lattice size on {R) for the Regge approach. The number A^o of vertices 
increases from 3^ to 8^. The well-defined phase survives an extension of the triangulated manifold 
which has almost no influence on the behavior of the computed expectation values {R). 
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